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Abstract. In this work, we use supratransmission and infratransmission in 
the mathematical modeling of the propagation of digital signals in weakly 
damped, discrete Josephson-j unction arrays, using energy-based detection cri- 
teria. Our results show an efficient and reliable transmission of binary infor- 
mation. 



1. Introduction 

The numerical study of wave transmission in long Josephson structures subject 
to harmonic driving was initiated in the middle 1980's by Olsen and Samuelsen 
[T], and the investigation was extended to coupled arrays of short superconducting 
tunnel junctions [2 . The study of the nonlinear bistability in Josephson junctions 
was the central topic of research in these works, a study that was continued later 
on via perturbation analysis with partially satisfactory results ^ until, finally, the 
whole analytical apparatus for the undamped, continuous-limit case was recently 
unveiled [4]. Nonetheless, the investigation in the area has proved to be rather rich 
and interesting [3 [6l [H E El [10] . 

In this paper, we develop an application of nonlinear supratransmission and 
infratransmission to the propagation of binary information in weakly damped, dis- 
crete, finite Josephson-junction arrays, by modulating the amplitudes of the driving 
signal at one end. The first section of this letter introduces the mathematical model 
under study and the energy expressions associated with the problem, while the next 
section briefly presents the finite-difference schemes employed to approximate the 
model. A computational study to predict the occurrence of nonlinear supratrans- 
mission under the presence of external damping is carried out in the next section. 
Next, we study briefly some characteristics of the kinks emitted by the boundary, 
and the following section describes the technique employed in order to transmit 
binary information in discrete Josephson junction arrays, presenting a simulation 
of transmission. Finally, we provide a section of concluding remarks. 



2. Mathematical model 

Throughout this paper, we assume that 7 and c are nonnegative real numbers. 
Likewise, we consider a system {un)n=i of TV Josephson junctions coupled through 
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(1) subject to : 



superconducting wires, satisfying the discrete, boundary-initial-value problem 

Un - C^A'^Un + 7„U„ -|- sin M„ = /Lt (1 < n < iV), 

' u„(0) =0, l<n<N, 

^(0) = 0, l<n<7V, 

(j)(t) 

Uo-Ui = —^, te(0,+oo), 

un+1-un = 0, te(0, +oo), 

where 7n = 7 for every n < N, and 7jv = 7 + The number R is called 

the output reading resistance of the system, and it is related to the output current 
intensity / through Ohm's law: / = un/R. The number c will be called the cou- 
pling coefficient, and 7 evidently plays the role of an external damping coefficient. 
The constant /i will be called the Josephson current of the system by abusing 
the standard terminology, and its inception in our model corresponds to a mere 
mathematical generalization of the model describing discrete arrays of Josephson 
junctions. Indeed, if /x is equal to zero, then <^(t) and I represent the input intensity 
function and the output current intensity, respectively, normalized to the Josephson 
critical current. 

Notationally, u and il denote the first- and the second-order derivatives of u with 
respect to time, respectively, and A^Un represents the spatial forward-difference 
Un+i — Un, while = A^Ax, where the backward- difference A^ = Un — Un-i is 
employed. The function </> may be assumed to be continuously differentiable over 
(0, -|-oo) in general. The particular case when </> is a sinusoidal function will be of 
special interest. More concretely, we will study the case when the driving is given 
by the function (j){t) = Asm{ilt). 

It is worth noticing that our mixed-value problem may be approximated via a 
continuous, Neumann boundary-value problem involving a perturbed sine-Gordon 
model for strong coupling. Also, notice that if 7 = /i = 0, then the Hamiltonian of 
the nth lattice site is given by 

(2) ^"^^^ - '"")^] + 1 - cos w„. 

After including the potential energy from the coupling between the first two sites, 
the total energy of the system becomes 

N 2 

(3) E = Y,Hn + j{ui-uof. 

n=l 

A simple integration of this formula over a finite interval of time gives the total 
energy injected into the system over such interval. 

The limiting case consisting of the discrete Josephson-junction array with an 
infinite number of junctions is of particular interest. In this situation, it is easy to 
verify that the rate of change of the energy of the system with respect to time is 
provided by the formula 

(4) ^=<^(t)^o-7E(^»)'- 
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3. Numerical scheme 

First finite-difFerence scheme. We consider the system of differential equations ([T|), 
and a regular partition = <o < < • • • < thi = T of the time interval [0,r] 
with time step equal to At. For each k = 0, 1, . . . ,M, let us represent the ap- 
proximate solution to our problem on the nth. lattice site at time tk by m^, and 
let us convey that Stu^ = wJi^"^ - '^n~^^ that S'^u^ = Mf;+^ - 2u^^ + u^^^ and that 
^'x'^n — ""n+i ~ + In order to posses discrete energy expressions that con- 

sistently approximate their continuous counterparts, the continuous-time equations 
will be approximated by the discrete expressions 

2.2 fc In , k y{ut^)-V{ui-^) ^ 



u„ 



for every n = 1, . . . , A^. Here, V{u) = 1 — cos u at any time. 

Our finite-difFerence scheme is consistent with our mixed-value problem, con- 
ditionally stable, and possesses energy-invariant properties. More precisely, if the 
energy of the system at the kth time step is computed using the expression 

N-l 2 

(6) Eu^Y. + Uul^' - "o+')K - 4)^ 

n=l ^ 

where the individual energy of the nth lattice site is computed by 

(7) i?^- = \ ( ""^^; i^{ui\\ - - u^) 

^ 2 

then the discrete rate of change of energy turns out to be a consistent approxima- 
tion for the corresponding instantaneous rate of change, while the approximations 
of the individual and the total energies are consistent estimates of their continu- 
ous counterparts. Particularly, the following property of consistency in the energy 
domain holds for every k — 1, . . . , M: 

At ^'^^''^2M-^^[2At 

n—l ^ 

Second finite-difFerence scheme. It is important to observe that the energy func- 
tional Hn defined in ([2]) is positive definite; however, there is no guarantee that 
the corresponding discrete approximation given by ([SJ is also positive definite. In 
order to save that shortcoming, we present now the finite-difFerence scheme 

where JuJ^ = ""n^^ +'"n^^- The local energy at the nth node and time tk is approx- 
imated then using the discrete expression 

2 2 



jjk ^ Ifu'^+^-u'^ 



(10) 2 V At 



i/(«f+i) + y«) 
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which is always nonnegative. Moreover, the total energy of the system is estimated 
through 

(11) 



This new scheme is clearly consistent with respect to problem ([IJ in both the 
solution and the energy domains. Moreover, the identity 



(12) 



Ek - Ek-i 
At 



Sun - Su'l\ StUQ 

2Ai 



7 



IV — 1 ^ 



■' \2At 
1 ^ 



holds for every k — 1, . . . , M. 

Computational remarks. In order to simulate a discrete array of Josephson junctions 
of infinite length, we consider a finite system ([1]) of relatively large length N, in 
which each 7„ includes both the effect of external damping described above and 
a simulation of an absorbing boundary slowly increasing in magnitude on the last 
N ~ No oscillators. Specifically, we let 



(13) 



In 



1 + tanh 



2n - No + N 
2^ 



+ 1 + ^^Nin), 



where 5jv('^) is equal to 1 if n = A'^, and zero otherwise. In practice, we will let 
K = 0.5, cr = 3, iV = 60, and Nq = 50. 

Also, in order to avoid the generation of shock waves in a system like ^ subject 
to harmonic driving (f)(t) ~ Asm{nt), we opt for increasing slowly and linearly the 
driving amplitude from to its actual value A, during a fixed period of time about 
the initial instant. After that, the driving amplitude will assume the constant value 
A. 

Both numerical methods have been employed to produce the results in this arti- 
cle. Our results have shown that both methods produce results that are in excellent 
agreement for small values of At, an observation that leads to conclude that the 
results in the following sections are intrinsic to the problem under study and not 
scheme-dependent . 

Computationally, the implementation of these finite-difference schemes is accom- 
plished via Newton's method, which in turn leads to implement Crout's technique 
with pivoting to solve tridiagonal linear systems. The stopping criterion employed 
in Newton's method was reached when the relative error of subsequent approxima- 
tions was less than 1 x 10"^, and the norm employed was the infinite norm in the 
space R^. 



4. Analysis of supratransmission 

The process of nonlinear supratransmission consists of a sudden increase in the 
amplitude of wave signals transmitted into a nonlinear chain by a harmonic dis- 
turbance at the end, irradiating at a frequency in the forbidden band gap. The 
existence of a nonlinear supratransmission threshold of the energy administered 
into a finite array of Josephson junctions satisfying ([T]) for a harmonic function 
(p and driving frequency in the forbidden band gap region < 1, has been es- 
tablished in the continuous-limit case [4] and numerically predicted for a discrete 
system. Using our numerical scheme, a prediction of the occurrence of nonlinear 
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supratransmission can be approximated for every driving frequency in the forbid- 
den band gap of the system of equations ([T]) , by estimating the value of the driving 
amplitude at which a drastic increase in the total energy of the system is detected. 
In the continuous-limit case, the driving amplitude As at which supratransmission 
first starts is related to il through the relation 

(14) As = 2c{l-fl^). 

We follow next the standard methodology employed in [TT] , and consider several 
values of the damping coefficient. The results of our numerical computations are 
summarized as graphs of minimal amplitude at which supratransmission starts ver- 
sus driving frequency, and they are presented in Figure[TJ The diagram representing 
the undamped scenario exhibits a good agreement with the continuous prediction 
for high frequencies, and damping is clearly shown to delay the appearance of the 
supratransmission threshold for frequencies ft > 0.35. 

On the other hand, the process of nonlinear infratransmission or lower trans- 
mission, as opposed to supratransmission, consists in a sudden decrease in the 
amplitude of wave signals in a chain harmonically driven at its end. The existence 
of a nonlinear infratransmission threshold under which there is a drastic decrease 
in the energy injected into the system was established in _4 in the continuous-limit 
case. As a matter of fact, the infratransmission threshold Ai is less than the supra- 
transmission threshold As, and between the amplitudes Ai and As lies a region of 
bistability. In order to approximate numerically the value of Ai for a given driving 
frequency in the forbidden band gap with associated supratransmission threshold 
As, we chose a driving function (j){t) — A(t) sm(D,t), with amplitude function defined 
via 

(15) Ait) = As{l - e-*/^^) + (Ao - - e"*/^^), 
where T2 ^ ti. 

A finite time period of relatively large length is chosen, and a finite, undamped 
system consisting of 60 Josephson junctions with no Josephson current, a coupling 
coefficient of 5, and a driving frequency of 0.9 is fixed, for which the nonlinear 
supratransmission threshold is predicted to be equal 1.945, according to Equation 
(fT4|) . Several driving amplitudes Aq are employed, and the total energy of the 
system is computed for each amplitude. The results of driving the system with an 
amplitude function given by (ITSI) with ri = 10 and T2 — 100 are shown in Figure 
[21 in which a drastic change in the behavior of the total energy with respect to 
^0 is detected around the value 1.41. This value is identified as the nonlinear 
infratransmission threshold of the system for fl — 0.9. 

As in the case of nonlinear supratransmission, the critical amplitude value at 
which infratransmission occurs also delays with the value of 7. More precisely, the 
larger the value of the damping coefficient the larger the critical infratransmission 
value. 

5. Signal propagation 

5.1. Character of solutions. Consider a system of 60 coupled oscillators de- 
scribed by ([1]) with coupling coefficient equal to 5, driven at a frequency of 0.9 in 
the forbidden band gap, and an amplitude of 2, just above the approximate supra- 
transmission threshold. Figure [3] shows the time evolution of a kink moving away 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Frequency 

Figure 1. Bifurcation diagram of occurrence of critical ampli- 
tude versus driving frequency for problem ([T]) with c = 5, for 
7 = (solid), 0.1 (dashed), 0.2 (dash-dotted), 0.3 (dotted). The 
continuous-limit prediction (jl4p of the critical threshold is repre- 
sented by a sequence of plus signs. 

from the driving boundary. The exact location of the kink can be better deter- 
mined by studying the time evolution of the local energies of the 60 sites at the 
corresponding times. 

With that purpose in mind, Figure U presents the evolution of the local energies 
of the six snapshots in the previous figures. The location of the kinks is accurately 
identified, at each individual time, as the absolute maximum of the local energies. 
Moreover, a constant phase velocity of approximately 1.44 sites per unit of time is 
observed. 

Finally, in Figure [5] we present a simulation of the evolution of the local energies 
at the same instants of time as in Figured The same parametric values have been 
used in both figures, except that in this last graph we employed a coefficient of 
damping equal to 0.001. In this situation, our numerical results (not presented 
here for the sake of briefness) show that the supratransmission threshold for this 
value of 7 is still smaller than the driving amplitude 2. Indeed, our results show that 
there is propagation of energy in the form of localized modes; however, dissipation 
of energy is present too, as evidenced by the fact that the local energy density of 
localized solutions is slightly smaller when 7 is not zero. 

In both of the cases presented here, the group velocity of the kink produced by 
the driving boundary is approximately the same: 1.44. However, the signal of the 
dissipative case seems to be a little delayed with respect to the case then damping is 
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1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 

Amplitude 

Figure 2. Graph of total energy vs. drivmg amphtude m an 
undamped, semi-infinite, discrete Josephson-j unction array, driven 
at a frequency of £7 = 0.9 via expression (|15|) during a period of 
time equal to 6000, with a coupling coefficient of 5. 



not present. This observation is in perfect agreement with the fact that the driving 
amplitude A is slowly and linearly increased from to 2, in order to avoid the 
creation of shock waves, so that the critical amplitude for the case 7 = is reached 
before the critical value for a nonzero 7. 

5.2. Simulation. Let 51 be a frequency in the forbidden band gap of our problem. 
Assume that Bi and Bs are nonnegative values that are just a bit smaller than the 
infratransmission and supratransmission thresholds, respectively, both associated 
to the frequency fl. We define the seed of the system as the function 

(16) Is (t) = i [{B, - Bi) sin(m) + B, + B^] , 

for every i > 0. 

Next, we define the period of signal generation P as an integer multiple of the 
driving period. In our case, P will be equal to 20 times the period of driving. 

Let a be a positive number with the property that Is{t) + a overcomes the 
nonlinear supratransmission threshold for some values of t. With these conventions 
in mind, a single binary message (61, 62, • ■ • , bi) consisting of I binary bits will be 
transmitted into system ([1]). In general, for every m = 1,2, ... ,1, we define the 
signal function Sm by 

ahm sin(rit), if mP < t < {m + 1)P, 
0, otherwise. 



(17) S„,it) = 
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Figure 3. Development of a kink solution of ([T]) for discrete array 
of 60 Josephson junctions, with 7 = /-t = and c = 5, driven at the 
end with by a frequency ft — 0.9 and an amplitude A — 2 slightly 
above the supratransmission threshold. The snapshots were taken 
at 6 times equally spaced between 90 and 115, employing a step 
size of 0.2 and an absorbing boundary in the last 10 junctions. 

The input intensity function will be defined then by 

; 

(18) (bit)^Is{t)+J2Sm{t). 

■m—l 

Consider a discrete Josephson junction arrays consisting of 8 junctions with 
coupling coefficient equal to 2. A driving frequency = 0.9 will be employed, for 
which the values B, = 0.23, Bs = 0.41, a = 0.05, 7 = 0.002 and R = 10 will 
be used. The binary sequence '10111001011011101001' is to be transmitted into 



BIT PROPAGATION IN JOSEPHSON JUNCTION ARRAYS 



9 













































































"^00000001300000000000000^1-,/-,^ 





20 30 

Site (n) 









































































°°Oor. n° 

"ooooooo" 




Ooooooooooooooooooooooorv . .-. 



Site (n) 









































































°°°ooo9ooooooooo 


oo*^ 







X 3 
2 
1 










































































ooooo'^ 









20 30 

Site (n) 















































































1 °°°°OoooononoooQG 



Site (n) 



LO 4 
?3 
2 
1 







































































°Oonnn,o.-,^r,n^'^on™iOOOOOOOO° 







20 30 

Site (n) 



20 30 40 

Site (n) 



Figure 4. Development of the local energies in a kink solution of 
([1]) for discrete array of 60 Josephson junctions, with 7 = /x = 
and c = 5, driven at the end with by a frequency J7 = 0.9 and an 
amplitude A — 2 slightly above the supratransmission threshold. 
The snapshots were taken at 6 times equally spaced between 90 
and 115, employing a step size of 0.2 and an absorbing boundary 
in the last 10 junctions. 



system ^ by modulating the input current intensity and read by the corresponding 
output current. The resulting outcome of our simulations is summarized in Figure 
[6] for a normalized time. 

Nonzero bits are clearly identified with output intensities of absolute values 
higher than the corresponding intensities associated with bits equal to zero. More 
accurately, nonzero bits are completely characterized by the fact that, at some time 
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Figure 5. Development of the local energies in a kink solution of 
([!]) for discrete array of 60 Josephson junctions, with 7 = 0.001, 
fj, = and c = 5, driven at the end with by a frequency fl = 
0.9 and an amplitude A ^ 2 slightly above the supratransmission 
threshold. The snapshots were taken at 6 times equally spaced 
between 90 and 115, employing a step size of 0.2 and an absorbing 
boundary in the last 10 junctions. 



in the corresponding period of reception, the value of the intensity of the output 
signal is higher than a cutoff limit (in this case equal to 0.01). 

Several other experiments have been carried out for different values of 7. Here, 
it is interesting to notice that our experiments have produced results which are 
qualitatively in agreement with that presented in Figure |6l Particularly, it is im- 
portant to mention that there is a well-defined difference between the bits and 1. 
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Figure 6. Graph of output current intensity vs. normalized time 
of system ([ij consisting of 8 coupled junctions with c = 2, as a 
response of the input current function ([TBI) for the binary message 
'10111001011011101001'. The parameters 7 = 0.002, n = 0.9 and 
i? = 10 are employed. 

The only difference is the amplitude of the oscillations of the output current, which 
decreases as 7 is increased, as expected. 

6. Conclusions and discussion 

In this letter, we have proved, using numerical computations, that it is possi- 
ble to transmit binary information in discrete Josephson junction arrays using the 
processes of nonlinear supratransmission and infratransmission. In the absence of 
dissipative effects, our model (which is based on the modulation of amplitudes of 
source signals with constant frequency) has shown to be highly reliable for suffi- 
ciently long periods of single-bit generation. Moreover, our computations show that 
the general picture does not change much when weak damping is present; indeed, 
the only difference obtained between the conservative and the dissipative cases lies 
in the amplitude of the oscillations of the output current intensity. In practice, 
the knowledge of the speed at which localized solutions move through a discrete 
array consisting of Josephson junctions in parallel and the relation between driving 
amplitude, supratransmission and infratransmission thresholds and damping, may 
be fruitful in the accurate design of binary signal transmitters, as well as in digital 
amplifiers and signal detectors. 

It is worth noticing the similarities and differences of our results with respect to 
the corresponding Dirichlet boundary-value problem. The results derived in this 



12 



J. E. MACfAS-DIAZ AND A. PURI 



paper in the supratransmission analysis for the external damping are qualitatively 
similar to those presented in [11]. On the other hand, the fundamental structures 
for transmission of wave signals in the Dirichlet scenario involved the propagation 
of moving breathers |12[ 113] : meanwhile, kinks and anti- kinks have been generated 
in the Neumann boundary- value problem. 
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